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summary' 

An analysis and calculations on the stability of the free laminar 
boundary layer between parallel streams were made for an incompress- 
ible fluid using the Tollmien-Schlichting theory of small disturbances 
Because the boundary conditions are at infinity, two solutions of the 
Orr-Sommerfeld stability equation need not be considered, and the 
remaining two solutions are exponential in character at the infinite 
boundaries. 

The solution of the stability equation is obtained in powers 
of (-i/aR), where R is Reynolds number and a is disturbance 
wave number. With an asymptotic solution as a start, the stability 
equation is numerically integrated. The eigenvalue problem between 
a, R, and the disturbance phase velocity c was thus explored by 
trial-and-error process. 

The calculations show that the flow is unstable except for very 
low Reynolds numbers . The regions of stability and instability in 
the a, R plane were checked by obtaining damped and amplified 
solutions on opposite sides of the neutrally stable solution. 


INTRODUCTION 

Some of the classical problems that have been of interest to 
many investigators in the field of the stability of parallel flows 
are the stability of Couette-type motion, Poiseuille-type motion, 
and b oundary- layer flows. Couette flow between rotating cylinders 
was successfully treated by G. I. Taylor (reference 1). The prob- 
lem of plane Couette flow, however, has not been decisively settled. 
Plane Poiseuille flow was treated by Heisenberg (reference 2), 
Pekeris (references 3 and 4), Goldstein, and Lin (references 5 and 6) 
and highly controversial results were obtained. Poiseuille motion 
in a circular pipe was treated by Sexl (references 7 and 8) who, 
after an incomplete investigation, concluded that the flow was 
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stable. The problem of boundary- layer stability was treated by 
Tollmien (references 9 to 11), Schlichting (references 12 to 15) 

Lin, and others, and an experimental verification of the results 
was obtained by Schubauer and Skramstad (reference 16). 

In all of the previous investigations, the problem of the 
stability of parallel flows was treated for a case in which there 
was at least one solid boundary. In the investigation discussed 
herein, a flow that has no solid boundaries is treated, the case 
of the stability of the laminar, free boundary layer between two 
parallel streams of semi-infinite extent in plane flow. In the 
analysis, the fluid medium is considered incompressible. After the 
steady-state flow configuration is determined, the stability of 
that flow configuration is investigated. The solution is carried 
out for only the case of flow in which one of the streams is con- 
sidered at rest. 

Most of the analysis for this report was done before June 1948 
in the form of a doctorate thesis in the Mechanical Engineering 
Department at Massachusetts Institute of Technology. Acknowledg- 
ment is made to Professor W. E. Hawthorne of the Mechanical Engineer' 
ing Department at M. I. T. for his general supervision of the work, 
to Professor C. C. Lin of the Mathematics Department at M. I. T. 
for his valuable assistance and suggestions, to Dr. L. H. Thomas of 
the Watson Scientific Computing Laboratory at Columbia University 
for his assistance in setting up the problem for machine solution, 
and to the International Business Machines Corporation for donating 
the use of the Selective Sequence Electronic Calculator for the 
numerous calculations. 


SYMBOLS 

The following symbols are used in the theoretical development: 
x positional coordinate in direction of principal flow 

y positional coordinate perpendicular to principal flow 

t time coordinate 

u velocity component in x-direction 

v velocity component in y-dlrection 

TJT stream function 

£ vorticity 
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P 

V = n/p 

5 

U 1 

U 2 

B = 8U^/v 

n 

f(Tl) 

x*- = x/8 

y* = y/6 

u i 

t* = t T 

u* = u/u x 

v* = v/U x 
TJT* 

®(y*) 

«p(y*) 

<p»(y*) * 


mass density of fluid medium 

absolute viscosity of fluid medium 

kinematic viscosity of fluid medium 

characteristic measure of boundary- layer thickness 

free- stream velocity of one of the parallel streams 
(taken as characteristic velocity) 

free-6tream velocity of other parallel stream, 
u x > U 2 

Reynolds number of boundary- layer flow 
y/ Vvx/U-l 


function, defining form of boundary- layer stream func- 
tion for time- independent flov 

dimensionless positional coordinate in direction of 
principal flov 

dimensionless positional coordinate perpendicular to 
principal flov 


dimensionless time coordinate 

dimensionless velocity component in x*- direct ion 
dimensionless velocity component in y*-direction 

dimensionless stream function 
** - t(y») + < P ( y .)e la( * , - ct ‘ ) 

steady-state part of dimensionless stream function 

amplitude of disturbance part of dimensionless stream 
function 

dy* 
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.. . df 

*’(*)) - dn 

a wave number of disturbance (always real) 

a g eigenvalue of a when B -» co 

c dimensionless phase velocity of disturbance 

c eigenvalue of c when B — >eo 

s 

i •T-T 


Erfc z = 

k l> k 2* k 3 
Im 

exp 


2 

dt 

constants of integration 
imaginary part of 

base of Napierian logarithmic system e raised to 
power in parentheses following exp 



THEOEETICAL CONSIDERATIONS 
The Steady-State, Laminar, Boundary-Layer Flow 
between Parallel Streams 

The equations for the steady-state incompressible boundary- 
layer flow with no body forces and no pressure gradient over the 
flow field are 


du , du 

u 35 + 




du . dv _ 0 


( 1 ) 

( 2 ) 


From continuity considerations (equation (2)), a stream function ¥ 
can be so introduced that 


u 


v = 


dip 

5y 

d¥ 

35 
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If the form of ¥ be specified as 

¥ = f(rj) 

T) = y/ Vvz/ U i 

then 

u = U x f * (tj) 

V = (rif’-f) 

and equation (l) simplifies to 

f f .» + 2f»" * 0 (3) 


From the physical flow configuration (fig. 1), the boundary 
conditions are 


y— ► + a> u->U^ 

y — ► - 00 u -* u 2 

or, in terms of the variables of equation (3), 

T)— *■+ oo f •—> 1 


TJ— *■ - 00 



Because equation (3) is of the third order, a third bo unda ry condi- 
tion must be specified. The third boundary condition is arbitrarily 
selected as 



It is of interest to discuss the asymptotic behavior of the 
boundary-layer equation (equation (3)). 


When 


T ) — *CD 
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and 

i 

f(n)*-»T i 

Therefore, as a first approximation for large positive values of q, 
equation (3) can be written as 

qf • • + 2f ' " SO 


and it therefore follows that 

1 - f* 2 kjL Erfc ($) (4) 

Similarly, for large negative values of q. 



Equation (5) degenerates for Ug/U^ = 0. For that case, if it is 
assumed that when 


q — > - oo 


f ( q )->0 

and 


f (q) — ► - a 


then 

I 11 

f'(T|)— »k 3 e 

It is also of interest to discuss the position of the boundary 
between the fluid of both streams. This boundary must be a stream- 
line and also pass through the point (x = 0, y = 0). From the 
form of the stream function V, it can then be seen that V = 0 is 
the bo undar y streamline. In order for this condition to hold at 
any value of x, it is necessary that f(q) = 0 along the stream- 
line. The boundary streamline is therefore 
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The significance of the q Q line can also be demonstrated from 

momentum considerations. If, for a fixed value of x, the boundary- 
streamline is located at the momentum principle can be formu- 

lated as follows: 




pu(U x -u) dy + / pu(U 2 -u) dy = 0 


( 6 ) 


If f and T) are then substituted into equation (6), the follow- 
ing equation is obtained: 


,T >1 


Uq 




t'(l-t') dq + / f ' ^tf " f 7 ^ = 0 

J- 00 ^ 


where q^ is the boundary between the two streams. 
Then 

+ oo 


f (-) - (i - £) f(^) - 1 


f(-«) - / (f) 2 dq = 0 (7) 


However, 


From equation (3) 


Therefore, 


J ' (f * ) 2 dq = ff' - J ff’' dq 


ff" = -2f ' 


/ <f ' 


) dq = ff’ + 2 f' 


When the preceding relation is substituted into equation (7), 

f(«) Yl - ^f(qi) - ^ f(-oo) - f (oo) f'(oo) + f(-oo) f •(-<*>) 


2f ' ' (°°) + 2f"(-o} = 0 


( 8 ) 
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However, 

f"(oo) = f"(-oo) = 0 

Equation (8) therefore reduces to 

f(») [ww] - ft—) - f'(-) - (i - ^f(ii) - o 

From the asymptotic behavior of f (equations (4) and (5)) 

f(») [l-f(«)] =0 

f.(-») w- ~ = 0 

L U 1 J 

Therefore 

t(\) = 0 

and = q Q , as previously demonstrated. 

The Stability of the Laminar Boundary Layer 
between Parallel Streams 

The criterion used in the investigation of the stability of 
the laminar-flow configuration is the behavior with respect to time 
of a smal l periodic disturbance introduced into the flow field. If 
the disturbance is damped with respect to time, the flow configura- 
tion is said to be stable with respect to that disturbance. If the 
disturbance grows with respect to time, the flow is said to be 
unstable. If the disturbance remains unchanged with time> the flow 
is considered as neutrally stable. 

In reference 17, Squire demonstrates that two-dimensional dis- 
turbances are more destabilizing than three-dimensional disturbances 
for two-dimensional parallel flows. Only two-dimensional disturbances 
are therefore considered in this discussion. 

e 

The equation of vorticity for an incompressible viscous fluid 
in plane flow can be stated as 

+ + v &y = 


(9) 
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where 


t- 


dv 

bx 


du 

dy 


t t 

S t dt *x " ax 



At 


sl 

a * 2 a 7 2 


Now, let £, u, and v consist of steady-state (time -independent) 
and disturbance (time- dependent) parts. 


Then 


where 


i = tU,y) + £'(x,y,t) 
u = u(x,y) + u'(x,y,t) 

v = v(x,y) + v'(x,y,t) 


i' « l 


u* « u 


V* « V 


If the expressions for £, u, and v are substituted in 
equation (9), the steady-state terms cancel out as satisfying them- 
selves. Because the disturbance is considered very small compared 
with the steady-state flow, all terms nonlinear in the disturbance 
are neglected. Equation (9) then becomes 

S't + u £'x + «'*x + ^’ y + v %- = 

It is assumed that u' and v' are of the same order of 
magnitude and that and £' are of the same order of magni- 

tude. If it is now considered that the steady- state flow is of the 
boundary-layer variety and that the principal direction of flow is 
the x-direction, then 
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Equation (9) then reduces to 

+ u J' x + v'^y « VA f 

If a stream function of the form 

V* = $>(y*) + cp(y*) e ia ( x *' -c ‘ t *) 


( 10 ) 


( 11 ) 


is introduced into equation (10 ) after the dimensionless variables u*, 
v*, x*, and y* are substituted for the corresponding dimensional 

variables, the Orr-Sommerfeld equation is then obtained (refer- 
ences 18 and 19). 

($'-c)(cp”-a 2 q>) - ® ,M <P - - jjjj (cp ,V -2aV’ + a. 4 <P) (12) 


where 






f'(Tl) 


A more rigorous and detailed derivation of the Orr-Sommerfeld equa- 
tion may be found in reference 20. 


If the characteristic length 


6 is set equal to 



then 


y* = T| 


and 


$ 


III - 


d 3 $ 

dy * 3 


ii _ f 1 1 1 

dri 3 


The Orr-Sommerfeld equation for the laminar boundary layer therefore 
becomes , 

(f ' -c) (cp' 1 —o? cp) - f"'cp ■ - ^ (cp' v - 2a 2 cp" + a 4 ^) (13) 

For boundary conditions on equation (13) , it suffices to specify 
that cp is bounded throughout the flow field. 

Because equation (13) is of the fourth order, a set of four 
linearly independent solutions to the equation exists. 
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'P “ c i<Pl + c 2 <^> + c 2^3 + C 4 C P4 

The nature of the solutions cpj and cp 4 can he investigated by 
introducing the following transformation into equation (13): 

cp a exp 

Equation (13) then becomes 

(f'-c) [(g 2 +g’)-a 2 ]-:f" ' = - ^[g 4 + 6g 2 g' + 3g' 2 + 4gg" + 

g' " - 2a 2 (g 2 + g' ) - a^J 

If g is expanded as follows: 

1 _1 

g - (aE) 2 g Q + g 1 + (ctE) 2 gg + . . . 

the following set of equations is then obtained by equating like 
powers of aE. 

(f'-c) So 2 = - 1 So 4 

(f'-cKso' + SSoSi) - -fUg/g]^ + Sg^gp') 

The successive approximations can then be obtained without 
integration 


5 SO* 

" ' 2 So 
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Finally 


_ 5 

cp 3 S (f'-c) 4 exp J ^ 


li oK(f'-c) dy 


5 

4 


cp. = (f’-c) exp 


It should he noted that when 


-] 

Ur aK(f^-c) dyj 


T|— cp^ — > oo 


^ >+co q> 4 


os 


Because cpj and cp^ are each unbounded somewhere in the flow field, 
they cannot be considered in the solution of equation (13). 


qp and cp can be solved for in the following manner: 
X U 


Let 


/ i (*0 

* ’5 ('.a) * 


(14) 


If eqviation (14) is substituted into equation (13) and like orders 
of (- X) are equated, the following system of equations is obtained: 

(f--c)(<p (k) '^aV k b-f' , 'T (k) -<P (k - l),T -2a.V k '- 1 )'’ + aV*- 1 ) 


For k = 0, equation (15), after being rearranged, becomes 


<P 


( 0 ) 


t I 


9 

or + 


f"A (o) 

—c)* 


= 0 


(15) 


(16) 


.(0) 


The asymptotic behavior of cp ' can now be investigated. Where 

q — >±00 

fin 


f'-c 


-*0 
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Therefore, 

<p^°^( ±a> )— e aTJ + C 2 e“ a1 * 

Because the solution must remain hounded, 

Cq_ = 0 when T)— »+» 

Cg — 0 TJ — >-00 

In the region where q/ 0 ^ is exponential, qp^ is also expo- 
nential, as can he shown hy substitution into equation (15). It 

(it ) 

can thus he demonstrated that op i 8 also exponential in the 
aforementioned region. The proper boundary conditions on cp are 
therefore 


T ) — > - oo cp' — >a cp 

T] — >+oo cp 1 — >-aq> 


From the boundary conditions and the general form of the solu- 
tion, the secular relation formally stating the eigenvalue problem 
can be easily obtained. 


F 


(a,c,B) : 


C P^' (~°°) - acp^(-oo) cpg* (-oo)-acf^(-oo) 


•P]/ (+ 00 ) + aCf^(+») cpg' (+»)+cu^(+a>) 


=0 (17) 


From equation (17), it can be seen that only particular combina- 
tions of the parameters a., c, and B allow the boundary condi- 
tions on qp to be satisfied. 

For a rigorous treatment of the asymptotic expansions for qp. , 

and *P4> 8ee referenc © 21. 1 
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SOLUTION OF THE BOUNDARY LAYER AND OKR-SOMMERFELD 
EQUATIONS FOR THE LAMINAR BOUNDARY LAYER 
BETWEEN PARALLEL STREAMS 

The general plan of solution of the boundary layer and Orr- 
Sommerfeld equations is to carry the analytical methods only far 
enough to reduce the problem to a routine numerical solution. In 
the problem under consideration, it is necessary to investigate the 
asymptotic behavior of the solutions so that boundary conditions 
at finite values of the independent variable we known in terms of 
the infinite boundary conditions. After boundary conditions at 
finite points are obtained, the solution can be continued by ordi- 
nary numerical methods. 


Asymptotic Behavior of Steady-State 
Boundary-Layer Equation 

The following transformation of equation (3) ff" + 2f'"® 0 
is permissible without altering the form of the equation: 

q(x) = i f(i) 

x = atj + b 

If the discussion is limited to the case where Ug/U^ = 0, it is 

then possible to expand a solution for the transformed equation as 
follows : 

lx |x 

2 2 

q = Aq + A^ e + A 2 e x + A 3 e +... 

The foregoing expansion is convergent for negative values of x. 
If a numerical solution is started from inside the area of conver- 
gence of the series solution, the asymptotic value of q' (») can be 
evaluated. Because the asymptotic value of f'(co) is 1, the scal- 
ing factor a is given by 


1 

/y/q'(®) 
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The scaling factor b can be fixed by designating 

f(il 0 ) = 0 

The asymptotic form of f is then given by 


2*n 


at) 


2^ 


^ = &q + a^e + a2© + a^© + • • • 


and, since f ' (starting) < c. 


f ,,t . 2 aT) . atj 2 aT) 

fT^ " V + V + V 


The various coefficients of the expansion can easily be evalu- 
ated by elementary methods. 


Asymptotic Behavior of Expanded 
Orr- Sommerf eld Equation 

If the asymptotic form of f' "/(f- c) is inserted into equa- 
tion (16), the following equation is obtained: 

(°)" ( 2 . ^ • aq \ (0) 

cp . - yc^+bje + b 2 e + b 3 e + . . J q, = 0 

and the form of follows: 


9 (0> . i + I 

(18) 

If cp^ is substituted into equation (15 ) to obtain the next 
order term cp , the following equation is obtained: 
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<P 


(i) 


t t 


•s® 1 ! 


- \a 


: +b, e" + boe aT1 + b,e + 


3 

o al l 


( 1 ) 


cp 


- l l e 


(oh- |a)ti 


+ l 2 e 


(cH-a)ri (ch- |®)il 


+ V' 2 


and 


(1) at) (oh- •sa) 1 ! (cH-a)T) o'" 1 ' 1 j. 

cp = e + h^e 2 + hg^e + h 31 e * + . 


(ch- |a)rj 


(19) 


From equations (18) and (19), the boundary conditions can he 
obtained for integrating the stability equations starting at a 
finite q. 

A tabulation of the values for the coefficients of the first 
few terms of the expansions is given in the appendix. 


Integration of Boundary-Layer Equation 
The boundary- layer equation is 

ff" + 2f'" = 0 


(3) 


or 


f'»' - - | ff” 


Differentiation once gives 


f ,v = - | (ff"» + ff") 


and twice. 


t 7 = - i (ff ,v + 2f'f'" + f" 2 ) 



Expansion in a Taylor's series gives 
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The values of cp^ * ' (il+w), cp^ (ti+w), Cp^ (n+w), and 

5cd ( 0) ,v / \ can now computed in a manner similar to the case 

( 0 ) 

of the boundary- layer equation and the integration of cp can 

(0) ,v 

he carried forth. As before, 6cf> can be evaluated in starting. 

Expansion of cp^ in a Taylor's series gives 

(j> (1) (i) + w) = ^ l) (t|)+wql l) '(il)+!|ll) (l) "(t|)+|- ' (n) + 


Now 5 2 qA')(tl) = Bcp^^q+w) - Bcp^^q) 
therefore 

5 2 cp^ 1 ^ (ll) = Cp^(t)+w) + Cp^^TJ-w) - 2cp^(tl) 


'.^C (l) "w + ^cp ( 1 >" , h) + ^cp (1 ) T, (,) + . . . 


and 


8 2 cp^ 1 ^"(Tj) = w 2 cp^ ,V (T|) + ^2 C P^ 1 ^ V, ( T1 ) + • • • 

t 

Therefore, 


cp (x) (n) - f^cp (1), ’(n) 




The quantity cp^ can be integrated as indicated in a step-by- 


step manner* From the asymptotic form of cp^ 1 ^, efcp^ 

12 

can be evaluated in starting. ' 

When and cp^ have been integrated to a large posi- 

tive value of q, the boundary condition at r\ = +oo can be imposed 
to evaluate B: 


,<i> ?! JV 


t t 


■p^ (0> + (sj) «* l> 

cp' (+0tJ + oup(+oo) = o 
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Therefore, 



Although the foregoing method of integration was used in the 
solution of the problem, it was not the most efficient scheme. 
Originally, the problem was programmed for integration along a 
rectangular path through the complex plane; complex integration 
was necessary because of the singularity at the point where f ' = c 
in equation (15). It later became apparent that the rectangular 
path was disadvantageous and that a better choice was a path start- 
ing in the third quadrant of the complex plane and travelling in a 
straight line to a point on the positive, real axis (fig. 2). Most 
of the original program could be used for the straight-line path. 

In order to expedite the solution, the equations were therefore 
integrated along the new path in the manner previously indicated. 
The interval of integration w was taken as 0.2 + 0.05i. 

A more advantageous scheme for the straight-line path of 
integration of equation (15) is included herewith for completeness . 


Method of Solution of Eigenvalue Problem 

If sets of eigenvalues of a, c, and E are to have any 
physical significance, E must be real; a is taken real and c 
may be complex. The process then to be - followed is: 

1. For a fixed value of a and a set of values of c, inte- 
grate the equations and obtain the corresponding values of E. 

These values of E are usually complex. 

2 . For the value of E that . is real , the corresponding c 
is the desired eigenvalue. 

3. Eepeat the process until the secular relation (equation (17)) 
is explored. 
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RESULTS AND DISCUSSION 

The secular relation (equation (17)) for the case of the "free" 
boundary layer was explored for damped, neutrally stable, and ampli- 
fied disturbances. The stable disturbance corresponds to Im c < 0, 
the neutrally stable disturbance corresponds to Im c = 0, and the 
unstable disturbance corresponds to Im c > 0 (equation (ll)). The 
disturbance equations were solved for values of Im c = 0, Im c = 
0.05, Im c = 0.10, and Im c = -0.05 (table I). The results were 
then interpolated for Ira E = 0 to give the neutral stability curve 
and curves of equivalent degrees of damping and amplification, 

(table II and fig. 3). 

The curves obtained were continuous, which indicates a low 
effect of rounding errors in the numerical techniques employed. 

Some integrations were redone using an interval of half that used 
throughout the problem. The results of both integrations were the 
same to the fifth significant figure and indicated the low trunca- 
tion error. Because the parameter of expansion of the eigenfunc- 
tion cp (equation (14)) was (-i/aR) and only the first two terms 
of that expansion were used, the determination of the eigenvalues 
of a, c) and R was inaccurate at low values of aR. For the 
same reason, it was impossible to obtain the lower branch of the 
curve of neutral stability. For purposes of comparison with the 
free boundary layer, Lin's curve of neutral stability for Blasius- 
type flow is included in figure 3. (See references 5 and 6.) 

A comparison of the stability characteristics of the free 
boundary layer with Blasius-type flow reveals certain basic differ- 
ences. First, the free boundary layer is unstable for some dis- 
turbances at an infinite Reynolds number (tantamount to "inviscid" 
flow), whereas the Blasius-type flow is completely stable in that 
range. This type of instability is caused by the presence of a 
point of inflection in the velocity profile. Some implications of 
the presence of a point of inflection in the velocity profile were 
first investigated by Rayleigh (reference 22). Another salient 
difference between the two types of flow is the presence of a highly 
oscillatory solution (<P 3 ) in the general solution for the dis- 
turbance in Blasius-type flow, whereas no such solution exists for 
the free boundary layer. This difference arises from the doubly- 
inf inite boundary conditions for the free boundary layer as opposed 
to the singly- infinite boundary conditions for Blasius flow. The 
foregoing is apparent from the asymptotic forms of the solutions 
<P 3 and Cp 4 . As shown by the form of the solutions, viscosity has a 

second-order effect on ^ and whereas it has a first-order 
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effect on ‘Pj and <P 4 . The inclusion of <P 3 in the general solu- 
tion of the disturbance for Blaslus flow therefore indicates a 
greater effect of viscosity on the disturbance than for the case 
of the free boundary layer. As could be inferred from the fore- 
going, therefore, the effect of viscosity on the stability char- 
acteristics of the free boundary layer is apparent only at very low 
Beynolds number, whereas the Blasius flow stability characteristics 
are much more affected even at higher Beynolds numbers. 

The inaccuracies due to small values of the parameter aB can 
be avoided by direct integration of the Orr-Sommerfeld equation for 
those cases. The asymptotic solution of the entire disturbance 
function could easily be developed as was done for the expanded dis- 
turbance function, and the numerical technique of integration after 
starting from the asymptotic solution would correspond to equa- 
tion (2l). These further solutions should be performed as soon as 
more high-speed computing-machine service can be obtained. 


CONCLUSIONS 

It is concluded that the laminar boundary layer between parallel 
streams is an unstable-flow configuration except at low Beynolds 
numbers. The method of calculation of stability characteristics is 
successful for small absolute values of the parameter (-i/aB). 

In comparison with Blasius-type flow against a flat plate, 
instability occurs at much lower Beynolds numbers for the free * 
boundary layer than for the boundary layer against a flat plate 
with no pressure gradient. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 21, 1949. 
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APPENDIX - COEFFICIENTS FOE ASYMPTOTIC SOLUTIONS OF cp 0 AND q?i 

a = 1.23849316 
a Q = -1.23849316 
a 1 = 1.23849316 
ag = -.30962329 
a 3 = .0860064? 





**20 


b l h 10 + b 2 

„ 2 
2aa + a 


b l h 20 + b 2 h 10 + b 3 
b 3o = ; — 

3 c 

3aa + 2 a 
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Kg = (2a) 2 (a + | a^ 2 h 20 

% = ^ 3a ) 2 ( a + f a ) 2 h 30 
2i = <3 q % 


Z 2 - do Kg + d i Ki 

2 3 = do K 3 + d i Kg + dg Ki 


J 11 = 


b l + *1 


a 

aa + o 


*21 = 
^31 = 


b l h ll + b 2 + l 2 


2aa + a 2 

b l h 21 + b 2 h ll + b 3 + z 3 
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TABLE I - EIGENVALUES OF a, c, .AND B 


[inviscid solution (B— »»); c g = 0.58727198; 
and a g = 0.39537 .J 


a 

c 

E 

a 

c 

E 



Im c, 

. “T 

0 

T 1 

1 

1 


0.390 

0.584000 

.390 

.584550 

.390 

.584581 

.390 

.584600 

.390 

.585000 

.390 

.588000 

.390 

.600000 

.385 

.581500 

.385 

.582190 

.385 

.582204 

.385 

.582500 

.380 

.579960 

.380 

.580000 

.380 

.583000 

.375 

.575000 

.375 

.577700 

.375 

.577850 

.375 

.578000 

.375 

.580000 

.365 

.573600 

.365 

.573900 

.365 

.574050 

.365 

.574100 

.355 

.567300 

.355 

.570000 


292.54 + 

34.121 

305.89 + 

1.961 

307.06 - 

.011 

307.27 - 

1.251 

312.08 - 

28.001 

210.66 - 

190.291 

3.96 - 

88.131 

156.95 + 

. 11.591 

160.64 + 

.251 

160.76 + 

.011 

161.85 - 

5.171 

109.78 - 

.011 

109.83 - 

.321 

109.18 - 

25.381 

79.73 + 

11.711 

83.69 + 

.681 

83.81 + 

.011 

83.94 - 

.671 

84.40 - 

10.051 

57.40 + 

.981 

57.49 + 

.381 

57.52 + 

.071 

57.55 + 

.011 

43.68 + 

4.081 

44.18 + 

1.131 


0.355 

0.571000 

.345 

.568600 

.335 

.566000 

.335 

.567000 

.335 

.567300 

.335 

.568000 

.335 

.570000 

.325 

.565000 

.325 

.567000 

.315 

.550000 

.315 

.555000 

.315 

.563000 

.315 

.568000 

.315 

.567500 

.305 

.570000 

.305 

.570900 

.305 

.571000 

.295 

.576200 

.295 

.576400 

.295 

.576500 

.295 

.577100 

.285 

.585000 

.285 

.586000 

.285 

.586300 

.285 

. .586500 


44.25 

- 

0.011 

36.19 

+ 

.04i 

30.83 

+ 

.621 

30.76 

+ 

.121 

30.74 

+ 

.02i 

30.68 

- 

.371 

30.45 

- 

1.351 

26.92 

+ 

.681 

26.75 

- 

.021 

24.35 

+ 

4.661 

24.40 

+ 

3.391 

24.09 

+ 

1.291 

23.63 

+ 

0 

i 

22.68 

- 

1.691 

21.12 

+ 

.16i 

21.00 

+ 

0 

i 

21.00 

- 

.021 

18.67 

+ 

.021 

18.67 

+ 

0 

i 

18.63 

- 

.021 

18.56 

- 

,09i 

16.38 

+ 

.101 

16.26 

+ 

.021 

16.25 

+ 

0 

1 

16.21 

- 

.021 
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TABLE I - EIGENVALUES OF a, c, AND R - Concluded 


a 

c 

R 


Im 

o, 0.05 



0.315 

0.562500 + 

0.05i 

72.57 

+ 

0.95i 

.315 

.562800 + 

.051 

72.68 

- 

• 07i 

.235 

.540000 + 

.051 

31.19 

+ 

8.001 

.285 

.555000 + 

.051 

32.55 

+ 

.251 

.285 

.570000 + 

.05i 

28.87 

- 

7.54i 

.285 

.580000 + 

- .05i 

24.47 

- 

10.721 

.265 

.548000 + 

.05i 

28.89 

+ 

2.461 

.265 

.550000 + 

.05 i 

24.79 

+ 

1.901 

.265 

.560000 + 

.051 

23.74 

- 

.791 

.245 

.560000 + 

.051 

19.47 

+ 

1.111 

.245 

.570000 + 

.051 

18.26 

- 

.151 

Im c, 0.10 

0.265 

0.550000 + 

0.101 

79.17 

+ 

7.141 

.265 

.551700 + 

.101 

81.12 

+ 

.12i 

.265 

.552000 + 

.101 

58.03 

- 

1.291 

.245 

.545000 + 

.10i 

39.83 

+ 

.881 

.245 

.548000 + 

.10i 

39.72 

- 

1.881 

.245 

.550000 + 

.101 

39.45 

- 

3.731 

.225 

.542000 + 

.101 

27.74 

+ 

1.481 

.225 

.544000 + 

.101 

27.58 

+ 

,75i 

.225 

.546000 + 

.101 

27.32 

+ 

.031 

.205 

.550000 + 

.10i 

21.06 

+ 

1.091 

.205 

.554000 + 

.101 

20.58 

+ 

.451 

.205 

.557000 + 

.101 

20.18 

+ 

.Oli 


Im c 

-0.05 



0.400 

0.500000 - 

0.051 

16.25 

+ 

14.251 

.400 

.550000 - 

.051 

26.91 

+ 

13.761 

.400 

.580000 - 

.05i 

34.34 

+ 

1.641 

.400 

.590000 - 

.05i 

33.20. 

- 

4.561 

.400 

.600000 - 

.051 

29.59 

- 

9.721 

.420 

.580000 - 

.051 

46.72 

+ 

10.091 . 

.420 

.590000 - 

.051 

'49.92 


3.101 
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